Efficient decoupling and recoupling in solid state NMR for quantum computation 
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A scheme for decoupling and selectively recoupling large networks of dipolar-coupled spins is 
proposed. The scheme relies on a combination of broadband, decoupling pulse sequences applied 
to all the nuclear spins with a band-selective pulse sequence for single spin rotations or recoupling. 
The evolution-time overhead required for selective coupling is independent of the number of spins, 
subject to time-scale constraints, for which we discuss the feasibility. This scheme may improve the 
scalability of solid-state-NMR quantum computing architectures. 
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I. INTRODUCTION 

Nuclear spins have been proposed as quantum bits 
(qubits) for quantum computation, due to their good iso- 
lation from the environment,-'^ With the aid of the estab- 
lished technology of nuclear magnetic resonance (NMR) 
and the use of "pseudo-pure" states^iiSiS, liquid NMR at 
room temperature has allowed the largest quantum com- 
puters to date.**. 

In liquid NMR, a solution is used in which each solute 
molecule has several nuclei, each nucleus having a differ- 
ent Larmor frequency in the presence of a magnetic field. 
This allows selective addressing of the nuclei through ap- 
plication of resonant radio-frequency (RF) fields at the 
corresponding frequencies. Each synthesized molecule in 
the aqueous solution undergoes the same logic opera- 
tions. 

Within each molecule, logic operations are performed 
by successive application of single-bit rotations (by RF 
pulse sequences) and multi-qubit conditional operations 
(by time-evolution in the presence of coupling among the 
qubits) . The design of such logic gates relies on the as- 
sumption that the coupling among the specific qubits can 
be set at will to on (coupled) for construction of a logic 
gate, and off (decoupled) for absence of a logic gate. A 
systematic scheme for this has been developed)^ but re- 
quires sequences of duration 0(n) as the number of qubits 
n increases, resulting in a clock-speed that decreases as 
the computer-size increasesi^ Because this slow-down is 
only linear in n, and because some quantum algorithms 
afford exponential speed-up, the scheme is efficient in 
principle. However, truly scalable operation will require 
quantum fault tolerancCf^ which requires logic gates to 
be sufficiently accurate and fast in comparison to deco- 
herence times to surpass the fault-tolerance threshold. 
Faster logic gates are therefore of crucial importance for 
future NMR quantum computers. 

As an alternative to a liquid, the use of dipolar-coupled 
nuclear spins in a crystalline solid'*'^ possesses certain 
merits for scalability, such as long coherence times^iS the 
ability to polarize the nuclear spins by optical pump- 



ing of electron spinsjii and more sensitive methods of 
nuclear-spin detectioni2iiS However, it introduces addi- 
tional complexity to the design of the pulse sequences 
needed for logic operations. In these proposed schemes, 
a magnetic field gradient is introduced to differentiate 
qubit ensembles, by analogy to the chemical shift in 
liquid-NMR. Unlike liquid-NMR, however, nuclei within 
a single qubit ensemble are dipole-coupled. The com- 
plexity of pulse sequence design arises from the additional 
burden of constantly decoupling these intra-ensemble nu- 
clear spins, which have the same Larmor frequency, while 
also decoupling and recoupling nuclear spins with distinct 
Larmor frequencies. 

In this article, we propose a scheme that addresses 
both of these issues: the linear slowing of logic opera- 
tions and the complication of added couplings in a crys- 
talline state. The required ingredients are that all spins 
must be coupled by strictly dipolar (or pseudo-dipolar) 
couplings in a large magnetic field, the field gradient 
or chemical shift distinguishing these qubits must be 
very large, and the couplings between different qubits 
must be much stronger than couplings within an ensem- 
ble. As previously discussed)^ this last condition may 
be achieved in a geometry where nuclei are placed in 
one-dimensional atomic chains, with a field gradient par- 
allel to the chains. Multiple, well-separated chains form 
the weakly coupled ensemble. Other geometries, such as 
well-separated ensembles of planar monolayers with two- 
dimensional gradientSfi^ are possible as well. 



II. DESCRIPTION OF THE SPIN SYSTEM 

In the presence of a static magnetic field along the z- 
axis, the Hamiltonian consists of the Zeeman energy of 
the nuclear spins and the coupling between nuclear spins 
by the dipole-dipole interaction. The Zeeman term is 
written 
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where 7 is the gyromagnetic ratio of the nuclear spin 
and is the magnetic field at the location of nuclear 
spin k. Due to the field gradient, the magnetic fields 
Bk and therefore the Larmor frequencies tOk's of the nu- 
clear spins differ. We denote the gradient-induced res- 
onant frequency separation between adjacent nuclei as 
LUk+i — oJk = 5uj. The average Larmor frequency, which 
we notate wq, is much faster than all other frequencies in 
the system. 

The dipole-dipole interaction between two nuclear 
spins k and I separated by the vector Vu is given 



Hr 



^0 ^^7^ 



47r 



' kl 



Ik li 



, (Ifc • Vkl){ll ■ Vkl) 



' kl 



(2) 



In the presence of a large applied field, the terms similar 
to I^If or I^I^ oscillate at frequencies near or 2a;o, re- 
spectively. These nonsecular terms contribute negligibly 
to the evolution in the timescale under consideration^ 
With this approximation, the dipole coupling for a pair 
of nuclear spins k and I takes the form 



When broadband pulses are applied, their time scale 
is l/ndoj <C l/(5w. Thus the dipole coupling between any 
pair of nuclear spins is well described by Eq. Q. For 
simplicity, perfect delta-pulses are assumed for broad- 
band pulses; each pulse can be applied instantaneously. 
In practice, the effect of finite pulse widths due to limited 
RF power in a real NMR experiments^ needs to be con- 
sidered and accommodated before the proposed scheme 
is implemented. We make a further assumption requir- 
ing greater scrutiny: we assume that selective pulses can 
be performed quickly in comparison to the free dipolar 
dynamics, thus allowing selective pulses that appear as 
(5-functions with respect to averages of the dipolar Hamil- 
tonian. This requires that 5uj 3> Dki- The applicability 
of these assumptions will be discussed in Sec. IIVI 



III. 



SCHEME FOR DECOUPLING AND 
RECOUPLING 

A. Broadband decoupling 
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and 9ki is the angle between Vki and the z-axis. The 
dipole coupling strength Dki between neighboring nuclei 
is typically several hundred hertz, but can be smaller 
depending on the geometry of the system. 

The nuclear spins are manipulated by an oscillating 
RF magnetic field Bi in the plane perpendicular to the 
z-axis. In the following, a "broadband" pulse (or "hard" 
pulse) with Rabi frequency 7-B1 = SIrf has a pulse du- 
ration ^ l/ripiF rnuch smaller than l/nSu. A "selec- 
tive" pulse (or "soft" pulse) with Rabi frequency wpjp 
has a pulse duration much longer than I/Slu in order 
to selectively rotate spins of a single Larmor frequency 
ujk- Nearest- neighbor spins receive an erroneous rota- 
tion of order sinc(7r(5aj/2cL'R,F); these erroneous rotations 
can be corrected by pulse shaping techniques discussed 
elsewhere ji4 

The couplings between nuclear spins and "broadband" 
and "selective" RF fields neglecting counter-rotating 
components'^ are described by the Hamiltonians 
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for broadband (B) and selective (S) pulses, where — 
± i/^. We denote the angular frequencies and phases 
as to and (j) for a broadband pulse and and 0/ for a 
selective pulse targeted on spin I. 



To establish our formalism for analyzing decoupling 
pulse sequences, we summarize the celebrated multiple- 
pulse sequence for decoupling, commonly called WHH-4 
(or WAHUHA)i&. 

We begin by considering the effect of a broadband 
pulse on the spins in the absence of selective pulses. We 
consider the usual "rotating reference frame," which ro- 
tates at the RF frequency to. In this frame the spins 
undergo a time evolution governed by the total Hamilto- 
nian Hz 



H^i , written in the form 
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where 5ujk = t^fe — ^- Since Q.txf ^ \5^k\ ^ Dki, the 
broadband pulse causes a rotation about the unit vector 
n = [cos 0, sin 0, 0] by angle = r^Rp^ a rotation we 
consider as instantaneous. The unitary operator for such 
a rotation is = exp(— «0J'^). The dipolar-decoupling 
sequence WIIII-4 consists of four such broadband pulses, 
all with — 7r/2. For this sequence we notate such a 
pulsed rotation with </) = as P^;, with </) = 7r/2 as Py, 
etc. 

To understand the evolution of the spins during a fast 
sequence of such pulses, we transform our Hamiltonian to 
a reference frame that follows the pulses. This reference 
frame is referred to as the "toggling frame." A Hamilto- 
nian H written in this frame will be denoted H{t). If 
there were no dipolar evolution or field gradient, then 
spins in the toggling frame would undergo no evolution 
whatsoever, even though the Pa pulses are causing rapid 
revolutions in the rotating frame. Toggling-frame dy- 
namics can be observed by "stroboscopic" observation, 
in which measurement with a particular phase reference 
is performed only once per cycle, when the toggling frame 
coincides with the rotating frame. 
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At the beginning of a cycle, the dipole couphng Hamil- 
tonian of Eq. Q in the togghng frame coincides with Hf^ ; 
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When we apply a Px pulse aX t — t and a Py pulse at 
t = 2r, then the toggling frame Hamiltonian becomes 
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and 



7Id(2t < t < 3r) = P^PyH^PyPx = H^- 
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Similarly, a Py pulse ait = At returns the toggling frame 
Hamiltonian to 

i/D(4r < t < 5t) = PxPyPyH^PyPyPx = i/g • (ll) 

Finally a Py pulse at i 5t brings the toggling frame 
Hamiltonian back to H^; 

Hb{5t < t < 6t) = PxPyPyPyH^PyPyPyPx = ■ 

(12) 

The decoupling pulse sequence WHH-4 shown in Fig.Q] 
is designed so that the system evolves in the toggling 
frame by for time 2t, for time 2t and H^, for 
time 2r in one cycle. The total cycle time tc is 6r. 

We see that this particular WHH-4 cycle of form 



It is important to realize that dipolar-decoupling se- 
quences have been heavily developed in NMR over the 
past 40 years. By compounding these sequences with 
super-cycles, they can be made robust against finite pulse 
width, pulse errors, and higher order terms in AHT. 
More modern sequences such as BR-24i2, and CORY-482tt 
have routinely proven their effectiveness in reducing the 
timescale for dipolar evolution by three or more orders 
of magnitude. 

Among such pulse sequences is MREV-16, a varia- 
tion of MREV-82ii2^, one cycle of which is expressed as 
lZ,Y,XjlZ,Y,XjlZ,Y,XjlZ,Y,Xj. The cycle time is 
tc = 24t. This sequence is more robust against pulse 
imperfections than WHH-4, but it also contains an addi- 
tional advantage. The simple sequence |Z, Y, XJ averages 
offset terms proportional to /| , and Iji to 
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while the compound sequence MREV-16 yields the aver- 
age operators 
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In the case of MREV-16, and lie in the plane per- 
pendicular to /| , while this is not the case for any WHH-4 
sequence. This property is critical for allowing selective 
RF rotations during decoupling, as we discuss in the fol- 
lowing section. 
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causes a spin operator /J in the rotating-frame Hamil- 
tonian to go to Ij, then to /J, and then back again, 
in the toggling-frame Hamiltonian. We therefore no- 
tate this cycle as {Z, Y, Xj following the notation of 
Mansfield—. Several equivalent WHH-4 sequences ex- 
ist. For example, if P^ is replaced by Pg to make the 
sequence t, P^^t, Py,2T, Py,T, P^-^t, then we would find 
going to — to ly , and then back again, and so we 
notate this cycle as |Z, 

In average Hamiltonian theory (AHT), the unitary evo- 
lution in the toggling frame under this changing Hamil- 
tonian is calculated by a time ordered exponential, ex- 
panded via the Magnus expansion--. The zero-order 
term is simply the average Hamiltonian over the period 
of the sequence 
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For any WHH-4 sequence, the average Hamiltonian for 
the dipole coupling is 

iiu = 7^(2ri/g -I- 2riyg + 2tH^) = 0. (14) 
or 

It may be shown that all odd-order dipolar terms of the 
average Hamiltonian vanish as well, due to the symmetry 
of the cycle»i^ 



B. Selective rotations 

a. Application of MREV-16 during selective pulses 
We now show how to produce selective spin rotations 
on desired spins while decoupling the dipole coupling 
between them by the MREV-16 broadband decoupling 
pulse sequence. The result of such a selective RF pulse is 
to affect the evolution of the desired spin ensemble, while 
leaving all other spins unaffected. 

Suppose we weakly irradiate the spins by selective 
pulses as in Eq. lO, while simultaneously applying the 
broadband decoupling sequence MREV-16. The RF 
Hamiltonian for selective pulses Eq. © is written as 



(17) 

in the rotating frame at the angular frequency oj of the 
broadband pulses, where 5ijj[ = — uj. Since fJ^p 3> 
Eq. H17|l varies slowly in time as compared to dura- 
tions and intervals of fast broadband pulses in the decou- 
pling sequence. Therefore, -ffpp is averaged by MREV-16 
as if Eq. H17() were time-independent; 
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Here we have made use of the fact that 1^ = ^ *^fe 
averaged to (3J+ + by MREV-16. 

We next show that this allows selective rotations on 
individual spins in the presence of the Zeeman term, av- 
eraged to 



(19) 



Since \5L0k\ S> |wrfI is assumed, Hz dominates the time- 
evolution of the spins and perturbs that evolution. 
It is convenient to use the interaction picture. 
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to find secular part of i?RF- Using the transformation 
of into -^^ ^j^g interaction 

picture, one finds that the secular part of Eq. (|2()|l is 
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where angular frequencies of the selective pulses lo'i are 
chosen to satisfy 5uj[ = SuJi/'i. This Hamiltonian de- 
scribes selective rotations of the target spins. The si- 
multaneous application of MREV-16 removes the dipole 
coupling for all spins simultaneously: Hd = 0. Nonsecu- 
lar terms involving a resonant spin (/) oscillate with fre- 
quency 2(5w;/3; nonsecular terms involving non-resonant 
spins (k 7^ /) oscillate with frequencies {5uJi±5LOk) /3. This 
is similar to the principles of selective rotations without 
broadband pulse sequences; errors on non-resonant spins 
are neglected if the pulse bandwidth is sufficiently small. 
The influence of MREV-16 on this consideration is that 
the bandwidth must be chosen 3 times smaller for the 
same level of error. 

b. Other decoupling sequences during selective pulses 
The choice of broadband sequence is not limited to 
MREV-16, but for selective rotations only certain se- 
quences will work. To discuss these constraints, let us 
generalize this result to an arbitrary pulse sequence. We 
still require the assumption that the sequence is very fast 
in comparison to fc;', so that Eq. (jSJ is averaged to a term 
of the form 



Hz + Hl^{t) = - ^ hSujkC -Ik-Y.^'^- 
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For example, Eq. (|15l) shows that the sequence {Z, Y, XJ 
yields C = (l,l,l)/3, I = (2,-l,0)/3, and r? = 
(0,1,— 2)/3. The problem is that ^ and t] are not or- 
thogonal to Therefore Eq. (|22|l contains terms of the 

form cos{uj[t-(j)i){C-^)C-Ik- and sm{uj[t-(t)i){C-T])C- h- 
In the interaction picture, these terms are diagonal in 



the eigenbasis of Hz; they do not represent transitions 
between energy eigenstates, but rather a time-dependent 
fluctuation of the energy eigenvalues. As a result, simple 
truncation of these terms leads to a poor approximation 
of the dynamics. It is very difficult to achieve a desired 
rotation controlling only 0; and Wpp. 

The situation with MREV-16 is much simpler, as 
C • ^ = C • ^ = aiid so these terms are not present. 
The resulting average interaction Hamiltonian is exactly 
the same as if there were no decoupling pulse sequence, 
except that resonant offset frequencies are scaled by the 
factor C = 1/3 and the Rabi frequency is scaled by the 
factor 1/2. Arbitrary rotations can be achieved through 



simple choice of (f>i and 



■^RF' 



just as in normal NMR 



techniques24. Other, shorter sequences satisfy the con- 
straint C ' ^ = C ■ ''7 = 0; for example, the sequence 
lZ,Y,XjlZ,Y,Xj has C = (0,0,l)/3, ^ = (2,-l,0)/6, 
and T] = (0, 1, 0)/3. The non-orthogonality of ^ and t] in 
this case, however, means that the relation between the 
RF phase 0; and the angle of rotation achieved by the 
selective pulse becomes more complicated. The general 
relation is 
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So, for example, with sequence |Z, y, X] |Z, F, X], we 
have 



17, 



where cfiQ = tan -'^(1/4). A selective 7r/2 rotation of spin 
I about X may be achieved with = Scoi/S, (j)i — (j)o, 
and uj'^pt = 27r\/T7/5. 

c. Simulation of fast decoupling sequences We have 
tested the principles of such selective pulses in the pres- 
ence of fast decoupling techniques with a computer sim- 
ulation. In this simulation, two dipolar-coupled spins, 
labelled j and k, begin in a random pure state and are 
subjected to MREV-16. Meanwhile, a selective tt pulse 
following the ideas above is applied at one of the reso- 
nant frequencies Suj'y The time dynamics of this process 
are simulated in two ways. In the first, the spin under- 
goes all the rotations of MREV-16 (treated as perfect, 
zero-duration rotations) while also evolving according to 
Eq. (|17|l . Here the evolution is found by direct time-step 
integration of the Schrodinger equation. In the second 
simulation, the spin state at time t is taken to be the 
expectation from the average Hamiltonian: 

m)) = exp{-tHzt/n) exp{-iH^^^''''t/n) |^(0)) . (24) 

The fidelity between the spin states calculated by these 
two methods is calculated for different parameters DjkT 
and (Sujj — 5uJk)T = SujjkT. The simulation is run several 
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times, averaging over random initial pure states. The 
resuhing average fideUty is shown in Fig. |21 We find 
that if T is short enough in comparison to Djk and Sujjk, 
the sequence works extremely well, and that errors due 
to a dipolar coupling which is too large or insuflicient 
selectivity are independent, at least at the two-spin level. 



C. Selective inversion for recoupling 

The combination of a broadband decoupling pulse se- 
quence and selective pulses allows independent single- 
qubit rotations, as discussed above. The coupling be- 
tween a pair of nuclear spins also must be restored in 
order to perform two-qubit operations. We will demon- 
strate how to selectively recover the coupling between 
a pair of nuclear spins without affecting other nuclear 
spins. 

We begin by applying rapid, broadband 7r-pulses 
that eliminate the offset term Hz- In what follows, 
these pulses are always applied on the fastest possible 
timescale. The average evolution under these 7r-pulses 
is therefore governed by the full dipolar term = ^^d 
[Eq. ©J. We then manipulate this dipolar term with 
selective 7r-pulses. These selective pulses may be accom- 
plished in the manner described in Sec. IIII B al or by 
the usual soft NMR pulses if the dipolar coupling is suf- 
ficiently small. The broadband 7r-pulses must be sus- 
pended during these pulses to allow selectivity. 

We now introduce a new "second-toggled" reference 
frame in which these selective rotations are treated as sin- 
gle, (5-function rotations. "Free evolution" in the second- 
toggled frame corresponds to the average evolution dur- 
ing rapid broadband tt pulses only and the dipolar Hamil- 
tonian, which is unaffected by those broadband 7r-pulses. 

Working in this frame, we selectively rotate nucleus k 
around the x-axis by tt (Xk). Then the second-toggled 
dipolar Hamiltonian is 

^ X.H^X, = Dki - Illf + Ipy). (25) 

We further manipulate the second-toggled frame by ap- 
plying broadband 7r/2 pulses that correspond to the usual 
WHH-4 sequences, except spacing them by a time T cho- 
sen to be much longer than the selective pulse widths 
but much shorter than the dipolar dynamics. Thus, as 
in Sec. IIII Al we have 

Hi{ 0+ < < < T-) = HI, 

Si{ T+ <t< 2T-) = P^HfP,, 

Si{2T+ < t < AT-) = Pr^PyHlPyP,, (26) 

Hi{AT+ < t < bT-) = P^HIP,, 

Hi{5T+ < t < 6T-) = HI, 

which averages to 
1 



Hi^— {2THI + 2TP^H^P, + 2TP^PyH^PyP,) 



We refer to this sequence as the Wk{X) sequence. 

Three more subcycles are required. In the Wk{Y) cy- 
cle, we selectively rotate nucleus k around the y-axis by 
TT (Yk). The dipole coupling in the second-toggled frame 
is then 

H^ ^ YkH^Yk = (-24!/f + I.^/f - ) . (28) 

Again taking this evolution through a WHH-4 cycle with 
pulse-interval T, the second-averaged coupling averages 
to 

H2 = -^ {2TH^ + 2TP^HiP, + 2TP^PyH^PyP,) 



--D,,{2ii^ + iliy). 



(29) 



In the Wk{Z) subcycle, we selectively rotate nucleus k 
around the z-axis by tt (Zk). 

The double-toggled dipole coupling is 

H| ^ D,, {2im + Illf + Illf) , (30) 

which averages to 

H:^ = -^ i2THI + 2TP^H^P, + 2TPsPyH^PyP,) 



6T 

\D,,{riif + iiiy 



(31) 



over a slow WHH-4 cycle. Finally, in the Wfc(— ) subcycle, 
we do not apply any selective pulses to nucleus k. In 
this case, our subcycle is just a fast WHH-4 where a 
Wk{X) cycle would have its selective Xk pulse, combined 
with a slow WHH-4 cycle. The resulting second-averaged 
Hamiltonian is nothing more than the dipolar coupling 
averaged through a WHH-4 cycle, which is of course 
zero. 

In this way, three kinds of effective dipole coupling 
Hamiltonians for a pair of nuclear spins are generated by 
selective inversion of one or two of the nuclear spins under 
a slow WHH-4 sequence. The important observation is 
that the sum of these three effective Hamiltonians is zero: 



Hi -f H2 + H^ 



0. 



(32) 



= --D,,{2IlIf + llll). 



(27) 



This fact will be used to decouple unwanted couplings 
when the coupling between a specific pair of nuclear spins 
is recovered. 

Suppose we apply the Wk{a) subcycle to nucleus k 
simultaneous with the Wk{P) subcycle to nucleus I (a, (3 
= — , X, Y, Z). The resulting second-averaged effective 
coupling Hamiltonian between the two spins will be Hi, 
H2, or H^ if a ^ /3, or if a = Table ^ summarizes 
all of the possibilities. 

Now we show how to recouple two spins k and I while 
the rest of the spins, represented by m in the follow- 
ing, are decoupled from both spins k and Couplings 
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TABLE I: Effective Hamiltonian generated by selective tt- 
pulses under a slow WHH-4 sequence. For example, when 
sequence W{X) (See Fig. ^ is applied to Ik and W{ — ) to 
the effective Hamiltonian over the cycle is Hi, Eq. 1271 . 
When the same sequence is applied to both spins, the dipole 
coupling is the same as the natural dipole coupling Hamilto- 
nian [Eq. iPJ] and averages to zero over a slow WHH-4 
cycle. 





Wi{-) 


Wi{X) 


Wi{Y) 


Wi{Z) 







Si 


S2 


Si 


Wk{X) 


Hi 





Hi 


S2 


Wk{Y) 


H2 


^3 





Hi 


Wk{Z) 


Hi 


H2 


Hi 






TABLE II: Timescales. It is assumed in the present scheme 
that each timescale is much faster than the one below it. Ex- 
pected orders of magnitude for these timescales are shown in 
the right column. 



Rate 


Description 


rad/s 


UJo 


Average Larmor frequency 


10^ 




Rabi frequency for broadband pulses 






Number of qubits in gradient times gradient- 


m^n 




induced frequency separation 




5lu 


Gradient-induced frequency separation 


10^ 


LORF 


Rabi frequency for selective pulses 


10^ 


y-1 

c 


Inverse of super- WHH cycle time 


10 


Dki 


Near-neighbor dipolar coupling constant 


1 


Aso 


Intra-ensemble dipolar coupling constant 


0.01 



between ensemble nuclei (that is, nuclei with the same 
frequency) are also removed. Apply Wk{Z), Wk{Y) and 
Wk{X) sequentially to Ik while applying Wm{—) three 
times to /,„, as shown in Fig. ^ The effective coupHng 
over this supercycle, consisting of three slow WHH-4 cy- 
cles with cycle time Tc = 6T, is averaged to 

^ (tc^s + TcH2 + T,Hi^ = 0. (33) 

The dipole coupling among m spins is decoupled as well. 

To spin Ii, apply Wi{Z), WiiX) and Wi{Y) in this 
order. The dipole coupling between Ii and Im is again 
averaged to zero. The coupling between Ik and /; is, 
however, 

(0 T.Hs + t.hA = ^Dkilk-h, (34) 

In this way, the two spins k and / are recoupled ac- 
cording to Eq. H34(l while the rest of nuclear spins remain 
decoupled from spins k and / and from each other. 

The present scheme works only in zeroth order in AHT. 
Dipolar terms of order n scale as t^D^^^ for other spins 
in the system with couplings as large as Dki, includ- 
ing spins distant from k and / which are supposed to 
be decoupled. To compensate for this, we must make 
Tc <C D'^j^ and perform quantum logic over many cy- 
cles. The resulting phase error for decoupled qubits over 
a logic gate lasting t ^ D^t/ scales as [rcDki)"'- In the 
perfect-pulse limit, all odd order terms can be eliminat- 
ing by simply symmetrizing the sequence, so that spin k 
sees 

Wk{Z) - Wk{Y) - Wk{X) - W*k{X) - W*k{Y) - W*k{Z). 
and spin I sees 

Wk{Z) - Wk{X) - Wk{Y) - W*k{Y) - Wt{X) - W*k{Z). 

In the sequence Wk{a), the phases and the orders of 
the selective pulses are re-arranged such that W^{a) is 
time symmetrical to Wk (a) . We refer to this sequence as 
"super- WHH." Its cycle time is Tc = 6rc = 36r. 



IV. FEASIBILITY OF THE PROPOSED 
SCHEME 

Although the present scheme does not suffer from the 
0{n) scaling of previous constructions^^ it can only work 
if the various timescales of the problem are well separated 
and high power pulses are available. In this section we 
discuss the feasibility of these assumptions. 

The approximations we used are summarized in Ta- 
ble ^ We now discuss the approximations made, the 
physical constraints behind them, the error created, and 
potential improvements. 



A. High DC field: luo > ^rf > nSuj 

Laboratory magnetic fields rarely exceed 10-20 T due 
to the finite critical currents in existing superconductors. 
Obviously, it is beneficial to have as large an external field 
as is available for this scheme, as the Larmor frequency 
sets the ultimate physical timescale for the computation 
rate. If some qubits see an external field comparable 
to r^RF, then nonsecular terms become important and 
Bloch-Siegert shifts must be considered. Phase errors 
due to Bloch-Siegert shifts scale as (fiRp/wo)^. 

B. High RF power: Qrp ^ nSuj 3> 5uiJ 

Although the present scheme does not slow down as 
qubits are added, maintaining the availability of broad- 
band pulses becomes more difficult as the number of 
qubits gets large. However, the amount of RF power 
available is limited only by such issues as heat dissipa- 
tion and arcing. With careful engineering and the use of 
microcoils,^'^ Rabi frequencies approaching the Larmor 
frequency are attainable. Phase errors due to insuffi- 
ciently broadband pulses scale for the most off-resonant 
nuclei as (nSoj / VLji-pY ■ 

Since f^RF is ultimately limited by the Larmor fre- 
quency ljJq, it would seem that the present scheme puts 
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an upper limit on n, at least for reasonable clock speeds. 
Further, if the RF power remains constant and the dipole 
coupling is weakened as n is increased, an 0{n) slow- 
down is again incurred. However, even if finite power is 
assumed, n need not be the total number of qubits in 
the system. One could imagine instead a field gradient 
of finite extent, in which only a subset n in a much larger 
set of qubits N is well-distinguished by the field gradient, 
while all N lie within a finite bandwidth accessible by the 
available RF power. The entire register can be accessed 
by moving the gradient, which could be done by switch- 
ing of different gradient coils (as in magnetic resonance 
imaging), by mechanical motion of ferromagnets (as in 
magnetic resonance force microscopy) , or by more exotic 
means involving manipulation of local electronic hyper- 
fine fields. By limiting a moveable gradient to a small 
subsection of the qubit register in this way, sufficiently 
broadband pulses can be applied with finite power; the 
qubits poorly distinguished by weaker gradients are still 
decoupled by the present scheme, ft is important to re- 
alize that such concerns will only be important in very 
large quantum computers with thousands of qubits. 



C. High field gradients or chemical shifts: 

Large gradients or chemical shifts are needed so that 
selective pulses do not become prohibitively long. Very 
large field gradients are attainable by a variety of means. 
Gradients as large as ^f T//im with sufficient field ho- 
mogeneity for sustaining measurable qubit ensembles can 
be made with ferromagnets,^'^ The qubit-qubit frequency 
separation Sui can be made higher than 10 kHz by suit- 
ably separating the qubits. Errors due to selective pulses 
affecting nearby nuclei scale as sinc(7r(5ti;/2wRF)- Com- 
pensation schemes for such errors are possible^i^ 

Although an ideal symmetric sequence has a vanish- 
ing first order term, a first order term does arrive in 
the presence of pulse errors and finite pulse width. For 
the simple sequences proposed here, the first order terms 
linear in the pulse width can only be neglected when 
the pulse widths of the selective pulses are much shorter 
than the cycle time. First order correction terms scale 
linearly with l/o^RpTc. It is important to note, how- 
ever, that compensation for finite pulse widths in mul- 
tiple pulse sequences can be achieved to high order by 
various means ^ 



D. Slow dipole couplings: ^ 3> D^i ^ Mso 

The timescale limitations discussed above limit the 
dipole coupling strengths which can be manipulated by 
this sequence. As discussed in Sec. lHlCl n*'' order resid- 
ual terms in the super- WHH sequence result in phase er- 
rors of order {DkiTc)". Higher-order sequences can alle- 
viate this restriction, but in the presence of non-idealities 



it is challenging to eliminate all second order terms. For 
the numbers used in Table ^ the dipole coupling must 
then be slower than I Hz, which is possible for reason- 
ably separated nuclei, especially if 7 is low. Unfortu- 
nately, this results in rather slow clock times. Although 
coherence times can be extremely long^ifi the exponential 
speedup of quantum algorithms can only compensate for 
this large, constant slow-down over classical computers 
for very large problem sizes. Previous proposals for NMR 
computers^ are no faster, as the clock speed for such com- 
puters with several hundred qubits would also fall slower 
than a hertz. We re-emphasize the important advantage 
to the present scheme: there is no adverse scaling with 
n with the super- WHH sequence. Therefore, if the clock 
cycle can be improved with higher magnetic fields, larger 
gradients, and error compensation schemes, the scheme 
will remain scalable for large quantum computers. 

The final issue not resolved by the present scheme has 
been discussed elsewhere.^ Copies of qubits in a single 
ensemble are decoupled by the compound WHH cycles 
of the present scheme. However, when different qubits 
k and I are coupled, all of the spins in ensemble k are 
coupled to all of the spins in ensemble I, leading to 
cross-couplings which ultimately couple ensemble mem- 
bers. This may only be compensated for in the present 
scheme and in existing proposals by spatially separating 
ensemble members in order to ensure that qubit couplings 
are much stronger than couplings within each ensemble. 
This limits available couplings to near-neighbor couplings 
onlyiS 



V. CONCLUSION 

We have seen that the judicious use of hard and selec- 
tive pulses allows rotations of single-spin ensembles and 
selective couplings between ensembles. Single-spin rota- 
tions are performed by application of a weak RF field 
in conjunction with a rapid broadband decoupling se- 
quence. Couplings between ensembles are accomplished 
by altering slower decoupling sequences with selective ro- 
tations, and constructing supercycles of such altered sub- 
sequences which couple only desired nuclei. The feasibil- 
ity of this approach depends on the selection of a physical 
system that exhibits a hierarchy of time scales. These 
time scales must be widely separated to minimize the 
phase error of hard pulses, ensure the selectivity of soft 
pulses, minimize the impact of finite pulse widths, abate 
higher-order contributions to the average Hamiltonian, 
and mitigate inter-ensemble cross-couplings. In solid- 
state quantum computer architectures where these time- 
scale constraints are satisfied, the present scheme resolves 
two previously unresolved issues, improving the tempo- 
ral scalability of such architectures: the linear slowdown 
of logic operations; and the complication of added intra- 
ensemble couplings. 
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FIG. 1: Decoupling pulse sequence WHH-4^ Pa and Pa stand for -k/2 and — 7r/2 pulses around the a-axis (q = 
The system whose Hamiltonian is H evolves in the toggling frame hy Z = H for time 2r, Y — PxHPx for time 2r and 
X — PxPyHPyPx for time 2r in one cycle = 6r. In the case of the dipole coupling Hamiltonian of Eq. Q, the zero-order 
term and all odd-order terms in the average Hamiltonian vanish over the period tc- 




FIG. 2: Fidelity versus rSuijk and rDjk- The average wave function overlap between dynamics calculated exactly and dynamics 
calculated with zeroth order average Hamiltonian theory. This simulation averages over 3 random initial pure states, and uses 
cj,, = 0.7, o-^p = Suj,k/WO. 
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FIG. 3: (a) Sequence Wk{X): Within each time-interval of a slow broadband WHH-4 sequence (consisting of 7r/2-pulses Px, 
Px, Py and Py), selective 7r-pulses around the a;-axis (represented by Xu) are inserted at the beginning and at the end of the 
interval. Between the two selective yr-pulses in an interval, fast broadband 7r-pulses are constantly applied at a speed faster 
than l/5uj in order to eliminate the Zeeman term. The selective tt pulses are produced by the combination of a broadband 
decoupling sequence and selective pulses. Sequences Wk{Y) and Wk{Z) are obtained by replacing 7r-pulses around the i-axis 
(Xk) by ones about the j/-axis (Yk) and the z-axis (Zk), respectively, (b) Sequence W{ — ): No selective 7r-pulses are applied. 
The fast broadband yr-pulses and a broadband decoupling sequence are synchronized with the sequences Wk{X), Wk{Y) or 
Wk(Z) that are applied to other spins. 
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FIG. 4: Super- WHH sequence to couple nuclear spins Ik and /; while decoupling the rest of the spins in a system Im from 
both. In the sequence W^^(a), the phase and the order of the broadband 7r-pulses are rearranged such that W^ia) has a time- 
reflection symmetry to Wk{a)- Nuclear spins with the same frequency are decoupled by this sequence as well. The restored 
dipole coupling between Ik and 7; is {8/9)Akilk-li- 



